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FOREWORD 


Iti the process of creating this packet of applied programs 
for calculating the navigational parameters of artificial Earth 
satellites (AES), the methods of use, whose possihilities and 
organization are described in Preprint [8] , there arose a library 
of subprograms which have a sufficiently independent nature. She 
present work is devoted to a description of these subprograms. 

Because of the space limitation of a preprint, the description of 
the subprogram library has been divided into two parts. Included 
in the present preprint are the subprograms for transforming 
coordinates and time, for determining the position of the Moon 
and Sun, and for calculating atmospheric density on the basis of 
various models of the atmosphere and disturbances specified by 
anomalies of the Earth's gravitational field. In the library of 
subprograms of the packet mentioned above, these subprograms have 
indexes A-B. 

Ihe second part of the subprogram library (F-l) contains 
subprograms for the formation of the right parts of a system of 
differential equations for the motion of AES and for its integra- 
tion by Adam's method, and subprograms for calculating the values 
of various functions from the parameters of the AES's motion. 

The description of the master program and auxilliary subpro- 
grams, which guarantee the organization of information input, as 
well as the calculation and printing or recording on magnetic tape 
of an arbitrary set of navigational paramaters (NP), makes up the 
content of an independent preprint. 

There is a variant of the subprogram library for execution with 
binary precision of the calculations enumerated above. 

When using the suggested library of subprograms or even one 
program of this library, one must keep in mind that all constants 
which are encountered in separate subprograms are located in a 
common region and values are assigned to them by a preliminary 
reference to the subprogram CONST (p. 2.1). Therefore, before turning 
to descriptions of actual subprograms, it is necessary to become 
familiar with points 1.2 and 2.1. 

The systems ‘of coordinates used and accepted when describing 
the subprogram library of designations are introduced in point 1.3. 

The principle of the organization of the library of subprograms 
is described in p. I.l. 

The remaining points contained in the description of actual 
subprograms can be used independently of one another. 

The author of subpxugrams VKMA and DENS (D02) is M.I. Ytyskovskiy , 
and of subprograms DEG2 , DEG3, and DEG5 (EOI) it is Ye. Ye. Ryazanova. 
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Ihe subprograms ADEN, AMBAR, GRAV, and WOCAI are taken from [ 6] 
and arn tested and modified for the Ye .Ye. Ryaaanova’a electronic 
computer (BSSM-6). She remaining subprograms are those of the 
author of this work. 

Ihe author wishes to thank Yo.A. Chistyakova for help editing 
the tests of the programs for publication and I.V. Zaytseva and 
V.V. Smirnova for their help in preparing the manuscript. 


A15XII.IARY SUBPROGRAMS POR GAIOULAIING IHE NAVIGAIIONAB 
PARAMETERS OF ARTIFICIAI EARTH SATEXBITES. FORTRAN IV. 

* V. I. Prokhorenko 

CHAPTER. I. THE ORGANIZATION OP THE SUBPROGRAM /5* 

LIBRARY AND ITS FEATURES 


1.1 Organ laat Ion of the SuBpro Pi ram Library 

At the base of the organ issational procedure of the subprogram 
library are the principles of organizing a subprogram library which 
are accepted at this time, and which are used, for example in the 
joint Institute of Nuclear Research in Dubna. 

To facilitate the review, the library is broken down into 
specific logical groups, each of which has its own index (A,B,0...). 
These indexes in a sense do not coincide with the indexes used in 
the joint Institute of Nuclear Research library but are used for 
the convenience of describing the library presented below. 

The chapter names conform to the names of respective groups. 
Subprograms in each group are numbered (for example, AOI, A02,...). 
Each individual subprogram is described in 0 separate point, and 
sometimes several subprograms which are linked to one another are 
described in one point. ’"The names of points coincide with the names 
of subprograms. In such manner the indexing can be considered an 
inventory of subprograms which are grouped according to their 
mean ing. 

Besides this, in each section of the description of the sub- 
program library is a list of subprograms by their names in alpha- 
betical order (together with indexes by which one can find a 
corresponding subprogram) . In the last section is a full list of 
subprograms given by names. 

When describing each individual subprogram the following format /6 
is used: 

1. Function 

2. Structure 

Subprograms, subprogram-function, and the packet 
of subprograms 

Identifier (identifiers) of the subprogram which 
is the input for the user. 

Internal inputs (subprograms inaccessible to the user) 
Peripheral subprograms used (access to other subprograms 
of the library)’ 

Peripheral devices (input and output devices) 
Common units (G0r#10N) 


^Numbers in the margin indicate pagination in the foreign text. 
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3. Access 

4. Input data 

5. Results 

6. Usage of the region COKKOU 

7. limitations 

8. Emergency outputs 

9. Method or algorithm 

10. References 

11. Text 

Information on all of the above points is not contained in 
each description, and numeration according to these points is 
not strictly adhered to. 

In order to save space in the description of structixre (p.2) 
replies of the type "Peripheral devices not in use" or "Access to 
internal subprograms unavailable" are omitted. 

1,2 C onstants, dimensional variable s 

The proposed subprogram library is a complex of subprograms 
which have been developed on the basis of several general principles. 

All constants, dimensional and non-dimensional, which are /j^ 
used in various subprograms, are taken out of the common region 
(COMMON), and values are assigned to these constants by accessing 
the subprogram 00N3I (for the majority of constants, see Cables 
2. 1-2. 4) and the subprogram GONGR (for coefficients of anomalies 
of the Earth's gravitational field, see Cable 2.5). 


All dimensional constants which are originally given in the 
system of units kg, m, and sec, can be subjected to multiplexing 
with the scale factors EM and ESEC, which are given as actual 
parameters of the subprograms CONSC and CONGR. 


Che problem is that for various AES it may be necessary 
to conduct the calculations in various systems of units; kg, m, 
sec; 1000000 m, 1000 sec, and so forth. Che system of units 
chosen for calculations can be fixed by two scale factors: 

EM, ESEC — the number of meters in the chosen unit of measuring 
distance and the number of seconds in the chosen unit of measuring 
time. In the case that the system of units is kg, m, and sec. — BM=1, 
ESEO=l . 


Dimensional reference data, such as T, s, y, z, v, v^, v v^, 

a, SB and so forth, should be translated into the system of units 
which are fixed by the scale factors EM, ESGC, for which the scale 
factors from corresponding units of OOMMON can be used (see 
Cable 2.3). 


In these subprograms for which the descriptions do not contain 
indication of the system of units in which the dimensional reference 
data should be fixed and resulting in dimensional results, it is 
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iiupllfd that it is a system o£ units fixed by the scale factors 
EM and ESE0» She current moment of time is given by the date and 
Moscow time S, figured from that particular date. She date can 
be given as the calendar date or 13 the RJD, the relative Julian 
date (nee p. 1.2). ^ 

She time S is mnasured in seconds (or in units determined 
by the scale factor ESEO). Only in the subprograms HMSSEO (BO 3 ) 
and SEOHMS (B04) is the time 'f always measured in seconds. 

In several of the subprograms of the library (usually the 
subprograms of other authors, for example DENS (DOE), ADEN (Do3)) 
the dimensional reference data should be given in definite units. 
Ihia is discussed in the descriptions of the corresponding sub- 
programs . 


1.3 Sy stems of Roordinatos. lime and Desifcnations 

1.7 The following system of coordinates is used. 
Greenwich relative rectangular coordinate system 0„,,_ 
with the rotating Earth: 


coordinated 


the center 0 coincides with the Earth's center; 

the axis 0- coincides with the rotational axis 
of the Earth and is oriented towards the North pole; 

the axis 0 is oriented towards the point of inter- 
section of the Earth's equator and the Greenwich 
meridian; 

the axis 0^ completes the system to the right 

Ihe absolute rectangular quadrate (equatorial, stellar) systems 
of coordinates OXYZ: 


the center 0 coincides with the center of the Earth; 

the axis OZ coincides with the Earth's rotational 
axis oriented towards the North pole; 

the axis OX is oriented towards the point of the 
vernal equinox (at the current moment); 

the axis OY completes the system to the right. 

Oscillating system of coordinates (elements of the orbit). The 
elements of this system of coordinates are: 

d 

-the semimajor axis of the orbit; 

6- -eccentricity; 

I -inclination (tie angle of incline of the plane 
of orbit and tha equatorial plane); 

Q -the longitude of the ascending vertex of the 
orbit (calculated along the equatorial arc 
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from the direction towards the point of spring 
Gounterclcokwise) } 

iu , 0 f perigee (angular distance from the 

ascending vertex to the perigee); 

-the time of passage through the perigee, 

Ibe position of the satellite in orbit is determined by the 
argument of latitude uu (angular distance of the ASS from th^ 
ascending orbital vertex). 

In descriptions of the subprogram library, the following 
concepts are connected with the estimation of time C 43 . 

-the Gregorian calendar (GU) — contemporary 
fhr''nol'',f^v 

-the Julian computation of time — the system 
of continuous count of days from the beginning 
of the Julian^period, year 4715 to the New era 
January 1, 12" according to the Gregorian calendar. 
-JD — the Julian date, the number of days 
which have passed since the beginning of the 
Julian period. 

-In many of the subprograms the relative Julian 
date (RJP) is used, the number of davs which 
have passed since 1900, January 0.12" of 
ephemerls time, 

RJDaJ3)-2 41 5020.0 

-Ihe stellar local time C 2 ] at the given 
meridian (s) is the time calculated from the 
moment of the upper culmination of the point 
of the vernal equinox to any other of its 
positions. Stellar time is numerically 
equal to the hour angle of the point of the 
vernal equinox. 

3. For several of the quantitites more frequently encountered /lO 
in the subprogram descriptions we will introduce constant designa- 
tions (identifiers of these quantitites). 

RJD-relative Julian date; 

1 ,1R- Moscow time , calculated from 
a certain date in seconds (or in other 
units fixed by the scale factor ESSO); 

50- stellar time on the Greenwich meridian 

at midnight in Greenwich on the corresponding 
date; 

51- stellar time on the Greenwich meridian at 
the moment of time 1 ; 
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YA-array 2C, t, g, V„, V„ 

(alJiclute systems pf crordinatio)}^ 
ya-aieray containing Xi y, z, v , Vy, v 
(Qreenwich coordinate 8ysta.ii)p ^ ^ 

Y*»array containing coordinates and 
constituents of a vector o.T ve’^ocity in an 
arbitrary system of ooordinates; 

XA-array containing only X, t, ’ii 
XS-array containing only x, y, fij 
X-array containing AKS coordinates in a 
random system of coordinates? 

A-array containing elements of orbit 

Cl f0f c f Si2 f to f cXt f 

SB-ballistic coeffioiont; 

1-atmogpherie density . 

When the sise of an array is mentionf'd in the descriptions, 
instead of the words ’’array Y is reserved for 6 real values," 
it will bo written "array Y,*." 

When the siso of an area reserved in the noMKON/B/3 blook 
is mentioned, the numin'r 3 indicates the quantity of real values 
for which block B is reserved. 



onm'M a 

C0NSUAK2S (INDKX A) 


/II 


E.l Basic noni^-.eantg (A01»»nQHST) 


1, Function. Ihe subprogram C0NS!U dispatches to the common 
area (CDMM) the values of dimeuBional and non-dimensional con- 
stants which are used in the system of subprograms for computing 
navigational information which characterises the positl^*^ ^ the 
AHS, and translates the dimensional constants into an glv n 
system of units. In Sables 2. 1-2. 4 a list is introdvv. I of cor- 
resnonding constants, their standard designations, the values and 
dimensions in the units (kg, m, and sec), for an assignment of 
the needed system of units the following parameters are used: 

BM-number of meters in the unit of measuring distance; 

ESEG-number of seconds in the unit of measuring time, 

for example, if for computations the chosen units are kg, 
ra, and sec, then BM»1, and BSBCwX. In order to carry out compu- 
tations in the systems of units, kg, 1000 km, 1000 sec then it is 
necessary to let BM«1,000,000, and SSBOrIOOO. 

2. Structure. Subprogram QONSI, 

Ihe Common Units 

/CA£D/i ,/CA00 /j,/CA00A/^,/:A22/^ ,/qor/^,/corl/^ , 
/CGRS/^ , /COM /^ , /CRE /, , /CRZ /J , /CAE/^ , /GAEL /, , /CCLZ /, , 
/COMZ/, , /COMZP/^ , /CSDAY /^ , /CTA ,\ , /CKDM/^^ , /CDSJS/, , 
fiPt/f, , /CDEOR/, , /CHRAD/^ , /CEZ/^ , /CE 3 /, , /Cf-i /, , /CE6/^ , A:C60,j, 
ATCSdOO/, . /3EM4 , /CEV/, , /CESB , /CELB/, , /CERO/, , /CHA/^ . 


5. Access; Gall COUSI (EM, ESEO) . 

4. Raw data; EM, ESEG. 

5. Results, use of the area COMMON. 

The values of constants in accordance with fables 2. 1-2. 4 are dis- 
patched to the units of Common enumerated in p. 2. 
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2.2 OonfficlRTits of Expanslnn of thfi TOarth's 

Gravliaiional Pleld In Spherical !!FuY3C'{:i~tY3S (AO^-DONGR) 

1. PuTiction: Sulsprogram OOHGR dispabches tha values of 

ooeffioieuts of expansion of the Karth'o gravitational field in 
spherical functions (5) to the common region (COMMON) and translates 
those coefficients into the given system of units. Values of the 
coefficients are given in table 2.5. I'hese coefficients are used 
only in subprograms for computing anomalies of the Barth’s gravi- 
tational field . 

E. Structure. Subprogram OONGR. 

COMMON units: /BnONGR/^^g. 

3. Access: OAIL OONGR (EM, ESEO). 


4. Haw Data: EM, ESEO /IS 

5. Results, use of the area COMMON. 

In the unit OOMMON/BOONGR/ANM (273), BN14 (?,73) , in accordance 
with tahle 2.5, the values o£ a are dispatched to array ANM, 
and the values ^ ^^^^are dispatched to array BNM. 
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Sa'ble 2,1 Astronomical units, gravitational characteristics, parameters of the Earth's 
ellipsoid, angular velocity of the Earth's rotation. 
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Table 2.2 Constants usod in the measurement of time. 


! 

• >Sm*< 

;Units ; 

;identi -■ 
ifiers 1 

Type 

Values 

Content 

S 



^ J 

5 

b 

I 

^CSDAY/, 



86400 

1 days in seconds 

z 

tCTb/ , 

r3 

REAL 

10800 

3 hours in seconds 

3 


nm\ 

INT& 

0 


January i 


(3 

ia 

GER 

31 


February 


, 


mm 

59 

Interval of time 

Mar :h 





90 

In days from Greer 

■April 





.120 

wich midnight , Jen 

May 





15 1 

i, of the current^ 

June 





I3I 

year to Greenwich 

July 





212 

midnight on the 

August 





243 

first day of the 
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December 
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number of days i 

n the yea 
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36525 

, ■ 1 

Julian centuries in epheme- 
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ris days 

J 

• 

Table 2. 

3 Auxilliary constants, scale factors 

\L±^ J 

a 
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b 
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3 

4 
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12 

13 
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PI 

PID2 
PI 2 
DEOk 
HEAD 
El 
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EA 
£6 
ce& 

€56Q0 

EM 
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2.4. Characteristics of the standard five-layer model 

of the Earth's atmosphere: /?/ - - layer boundries 

by altitude, coefficients 

( I ‘='i, 2 '’,‘: .v, 5 ) 


' Units 
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Here ANM (l)=ot 2 o =0, in ss far as the corresponding member of the 
resoultion is considered to be in the Earth's normal field (the 
value a is from table 2.1). 
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CHAjP1}’>3R 5 

SUBPROGRAMS POR iHE 5RANSP0RMATI0N 
OP TIMS AND OCORDINATBS 
(INDEX B) 


3.1 Transit Ion from GU, the Calendar Date . 
to RJD. Relative Julian Date (BOl - 
! QAtDAY ) 

a. BwawMMPmaM* ' 

1. Punctlnn; the program determines the RJD, the relative 
Julian date, the numUer of days which have passed since the mean 
Greenwich midday, January 0, 1900 until the mean Greenwich mid- 
night of the given calendar date. 

2. Structure. Subprogram; DATDAY. 

Common units :/GKDM/^ 5 ,/G32/^. 

3. Access; CADI DATDAY (DTK, RJD). 

4. Raw data; DTK-real variable containing the date GU which 
is given by the decimal fraction of the type; 0. DDMMGG, where DO 
is the date. MM is the month’s number and GG is the year’s number 
minus 1900 (the two last numbers of the year). 

3. Results; RJD-rslative Julian date. 


6, Use of the GOMMCN units. Constants are used from the 
units; /CKDM/^^ ,/GEZ/^ (No. 3 Table 2.2., No. 4 Table 2.3). 

7. Algorithm; RJD=365Ga + L + t^^ere 

GG is the number of years minus 1900, the number of days passed 
since January 0 of the current year until the given date, [ x] is 
the whole part of x. /20 


8. Text 

S;itlP0Ur r ‘iR Oi ^ DaV ( Ci t f C ' 'J 
0 5M »f)'|/C 52 /Er; 
cnii.MO-i/nOM/KOhf U) 

'•< = OAr <»[;?. 

Unrns.i 

/ = < <-.|jT I ) •F.n 

i/rl 1 #^ |M J , (T^ 

3.2 Transition from RJD (the relativ 
Calendar Date (B02-DAYDAT) 


n.)Tr,= k;3.M( NM ) *‘i C’-n 
K u - ( •! r. “ 1 > / A 
I r { < ) 2 » 1. . 2 

1 I r ; ‘ I ' • ~ ) 2 f T 
3 K'lTtisKOTr.- 1 

HRTu/; I 

i ’ 

3 Julian date) to GU, the 


1, Punction. Por the moment of time, given by RJD the relative 
Julian Date (see p. 1.3), the calendar date is determined. 
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The moment of time can also he given in the form of RJD and 
the time T in seconds (or in units given by the snale factor SSBO, * 
calculated from midnight of the RJD date (the interval of time 
T may contain any number of days). A result of the operation of 
subprogram DAXDAT, the interval which is a multiple for whole days 
is excluded fr-’^m the time T; TR is obtained, and the days excluded 
from T are added to the date RJD and the RJDR is obtained. The 
calendar date is determined for the RJDR date. 

2. Structure. Subprogram DAYDAT 

Common units :/0SDAY/^ ,/OKDM/^j ,/OOGO/^. 

5. Access; CALL DAYDAT (RJD, T, RJDR, TR, ND, NM, NG) 

4. Raw data; RJD-relative Julian date. T-time (in seconds), 
calculated from the RJD date. 

5. Results; RJDR-relative Julian date; TR-time in seconds, 
calculated from the RJDR date; KD, NM, NG-whole variables containing 
the corresponding number, month and year (minus 1900), corresponding 
to the RJDR date. 


6. Use of the area COMMON. Constants are used from the /21 
units CGMMOK/CSDAYA ,/OKDM/t*./OOG0/t (No. 1.5 of table 2.2, No. 8~ 
of table 2.5). ^ ^ 


7. Algorithm. let RJD be the reference Julian date, T the 
time in seconds calculated from midnight of the RJD date. 


PMR = PaD+[7ySD/^/]; TP.-T-[t/iSDAV] • 


where SDAY is the number of seconds in the days, [x] is the whole 
part of X. 


KV6= (RJDR-365)/1491 ; K'sRJDR-KVG; N = K/365, 

where KVG is the number of leap years which have passed sine 1900 
to the given moment. 

NG= I ^ ^ “ 

t N-1 , if 4 (KVG+1 )-k < I 0. 

The number of days, calculated from the beginning of the 
present year; 

KDTG=K-565NG + 1. 


The month and day of an unknown date are determined with a 
comparison of KDTG with the array KDM- (i=-.1.15), in which the i 
element contains the number of days which have passed since the 
beginning of the present year (not a leap year) to the beginning of 
the i-raonth. In the case of a leap year, when XDTG ^60 a correc- 
tion is made. 
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3, Text. 

s I- i t )-'( ',4 J-;-? ( ;;flv * T , i 'h . S5 . ‘r. nC ! 

L ■) I ( 

f. I i ■. *. J »•)/ • { 1 'n 

r. *! ''.ri I i »M w,*$ ') • y 

K’.-f'i ,) A / 

y.^n 

T .* s T - . f, ) A U 

'..iv"* ! r « H / ; ♦•''c 

(,)•<-%*)•>! .’.15 5/.46X 
; s^^.OM-r.v"; 
fJ.'JSK I i ) 

I -''js ! .V*',- .1 > O A-*'n 
: '•j '.''I > ii • 2 
iJ 5 = N'5 - , 


imtr*im»if ^ * •*«► »'*•• • i>« * 


2 Kpi GsK* I- J5*t'nn{ ; 

1 r n V 3 ) ^ I > I 6 

3 I r< KlTG -C^>U ) 6, 5 

4 un = f!9 
NM~ 2 
sn’To i3 

r, .(oTCls ;D'3-i 
f. D ) ' .’:•■ = i < 13 

IF< iVjTj-KOht »!H) » »t t 
’ cnKrjrjuE 
'o •IMSHU-I 

N') = KOTO“K 5.'t : f,M) 

KETUf^N 

EMi) 

5 QQ»vorslQTi or Hou rs. Minutes 
and Seconds into Seonuds 
XBP5-HMSSl5tr ^ 

1. Function; Originating from the time given in hours, 
minutes, and seconds, the time is determined in seconds. 

2. Structure. Subprogram HMSSBG 

Common units; /CB2/j^,/COG0/^./OC3600/3 . 

5. Access; CALI HMSS3C (HMS, T). 

4. Raw data; HMS-the real variable, containing the time 
given in hours, minutes and seconds in the following form; 

0. HHt'lMSSDS, where HH is the hours, MI-1 the minutes, and SSLS the 
seconds with fractions. 

5. Results; T-the time in seconds. 

6. Use of the COMMON units. Constants are used from the 
units GOMMON/C3S/3^,/GCGr./^,/OO3600/^ (Nos. 8 and 9 of Table 2.3). 


/22 


7. Text. 


■i l*i 

C 1 I 

0 1 V 
■J 0 

* , 

1 ' ; ' ^ . 

' <r i ,1 * 


)-r I", Mass'^r. ' H,i3 
to = 

to f/cc.to.i-j/csooa 
.•t'fi/r.m u/r.<s'i 


r. n T :t 5 
A.’tSs</«»r p, 

/ts«ns- , i 

ft'l.SffAMC-'? 

! F t A I ^ j , ^ I 

A.ts?" , 

: * ' ' n - * c ^ 0 » « . , n 5 5 3 g 

* * ' * • I 

S !•) 
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3.4 Oonvegflloti of Seconds Into Houre . 

Minutes. and leconds 

(B'b-t-SSaHMS)' 

1, function; Originating from the time given in seconds, the 
time is determined in hours, minutes, and seconds. 

2, Structure: Subprogram SSCHNS, 

Gommon units: /OGG0/2^,/GG 3600/ 

3. Access; OAII SB0HM3 (l, m, KM, «iSC). 

4. Raw data: I-time in seconds and fractions of seconds. 

5. Results; KH-hours (whole variable), KM-minutes (whole /23 
variable) , SllO-seconds and fractions of seconds (real variable) . 

6, Use of the area COMMON. Constants are used from the 
units GOMMON/OOG0/3^,/C3O3600/^ (Nos. 8 and 9 of Table 2.3). 


7. Text, 


r [ ,% !, ■: - SEv S 

STM , "•S' ", rfio 

. *.-*■ 'f.S '^C3/C?6 0; 

•: H 5 T / f 1 D ’< 

-t 5 !; = T„| 

i 

I ?• i / . i . . 1 . 

* I* ♦» I 



3.5 Stellar Time (S03-STTIME,, 3TT) 

1. Function: Subprogram STTIME determines the stellar time 

SO on the Greenwich meridian at Greenwich midnight of the given 
RJD date. 

Subprogram STT according to the known SO determines stellar 
time ST on the Greenwich meridian at any moment of time T on the 
RJD date. 


2. Structure; The independent subprogram STTIMF and the 
subprogram-function STT . 

Common units: /CSDJS/^,/OPl/3,/CCME/.^,/CT3Aj^,/BSO/5. 
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DescriptioTi of the subprogram S!0]’IJ4E. 
Access; GAIJi SflJIMS (RJD, SO), 


Raw data;RJD-relative Juliah date. 

Results; SO-stellar time ou the Greenwich meridiau in 
radians ( 50 Ss 2 ,11 ) . 


Gse of the units OGMKON, Constants are used from the units /24 

COMMON/CDS JSA , /CM/, (No. 4 of Table No. 1 of Table 2.5). ~ 

✓ 


Algorithm; 

so - e''38"45’ ,836 + 8640184*, 542 T + o‘,0929 T® + o‘,06II64Ai|), 
^iv|) s W.T7 123 Sin f 

.Si a 259°I0'59“,79 - I934°08' 31" ,23T + 7",48 t2 + 0" ,00801®, 
r ~/iJ£> /36525,/?JZ) 


4. Text, 

*?Ue'^0UTlNE STT ImE 1 OdVf 50 ) 

COMMON/COS JS/CSJ5 
COMMO^’/OPJ / Pf *P| !;2 , p J2 
TsOA^/OSjS 

S0m-SIN< I < .3®765lE-7*T^,3«2fi4c634£-<,-3j,757lt625)*T 
♦ . ♦<•,5236 0 I5i6) > 

W56?6,33#t 
N = .V / P I 2 
Wsrt-N*P ! 2 

SOs(6,7.'5 5F-6#T*.19sUF.-2J<'T*i.7399 3 5<J-.766385p-<,,sOW*W 

NsS0/pl2 
S0sS0-\'*Pl 2 
*1“ T'JRN 
END 

5. Description of the subprogram^f unction STT. 

Access; ST»STT(T). 

Raw data:T is the present moment in time in seconds, 
Moscow time. 

Results: ST-otellar time on the Greenwich meridian in 

radians. 

Use of the region COMMON; in the unit OOiWN/BSO/SO, TS, 

NS the fi^llowing values should be initially placed; 

SO-stellar time at Greenwich midnight on the date RJD, 
TS»0, SS«0, if T is calculated frnm the same RJD date. 
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In most eases ths RJDI (late and ths looking of time I can /25 
be Tes© -fchan or equal to RJ3. If RJDI < RJD, then IS and NS ^ — 

should contain (IS in seconds, and NS in days) th<^ interval of 
time between the RJDl and NJB dates. 

Constants are used from the units C0MM0N/00M2 /t ,/GI 3A (No. 14 
of Table 2.1 and No. 2, Table ?.,2). ^ ^ 

Algorithm. Stellar time SI is calculated according to the 
following approximate formulas 

5r=5o+«j(7'-r&-t08oo‘) 

Where is the absolute angular velocity of the Earth‘s 
rotation. 

SO-is the stellar time at Greenwich midnight on the RJD date, 

I- is the Moscow time, figured from the RJDI date, usually 

different from the RJD (RJDI < RJD) date, 

IS-86400® (RJD-RJDI). 

lext : 

STTfT) 

, TS , fis 

cnM-n i/co.'U/r/ij 

S 7 T S •) * r, I c - T > ) 

«STii5;’j 

K-|0 


5 • 6 Iransitlon from the Abi lute System of Coordinates 
to the Greenwich and the Rever se! 
from the G'reenwioh System t o tTii Absolute 
(b66-aGiga .agigag ) 


1. Function. According to the known values X,Y,Z,Vj^,Vy in 

the absolute system of coordinates, the subprogram AGIGA determines 
the values of x,y ,z,Vj^,Vy,v^ in the Greenwich relative system of 

coordinates; the transition back is also possible; x,y ,z,Vj^,Vy jV^ 

X,Y,Z,V ,V.,V„. The subprogram AGIGAG makes the same conversion 
X y z 

possible only with the coordinates: X,Y,Z ^ x,y,z and x,y,z’-^ 

X,Y,Z. 

2. Structure. Subprograms: AGIGA, AGIGAX 

Common units: /COME/^,/CT3/^. 




5. Access to AGIGA; CALL AGiaA (SO, T, Yl, 1,YR). /ZS 

Raw data; C-Moscow time (in seconds), 

SO-stellar time on the Greenwich meridian at Greenwich mid- 
night j Ylg-array of reference values of coordinates and consti- 
tuents of the velocity vector; I- indicator of the transition, 
lal during the transition from the absolute system of coordinates 
to the Greenwich, during the transition from the Greenwich 
system of coordinates to the absolute. 

Results: YRg-array of values of coordinates and constituents 

of the velocity vector in the resulting system of coordinates. 

4. Access to AGIGAC; OAll! AGIAC(S1D, X1,1,XR). 

Raw data; St-stellar time on the Greenwich meridian at 
the present moment of time; xl,~array of base coordinate 
values; 1-indicator of the coordinate transition (which 
both by its meaning and values conforms to the 1 parameter 
in subprogram AGIGA). 

Results; XR — array of coordinate values in the resulting 
system of coordinates. 

5. Use of the area G0r#I0N; in subprogram AGIGA constants 
are used from the units COMMCU/COMZ/t ,/0T5/t (No. 14, Sable 2.1, 

No. 2, Table 2.2). ^ 

6. Algorithm: a) transition from X,Y,Z,V ,V ,V to x, y,z, 

v^, Vy, Vg, is conducted according to the formula; ^ ^ 

a:=X cosj5 + Y sin p , A r 

y- -X sin (5 4 Y cosp , sinp+ Vy cos^-oy^jc , 

Z-Z, ^2 * ■ 

p .= 30 + - s'* ) , 


SO-stellar time at Greenwich midnight, 

V/3-angular velocity of the Earth's rotation, 

T — Moscow time /- 2 I 

b) transition back from x,y,z,v ,v ,v to X,Y,Z,V V ,V_ is 
conducted according to the formula; ^ ^ ^ x y z 

X cos^ -y sin |5 , = ^cos^ sinp» - , 

sin^ + </cos^, Vy sinp + x;^ cos^ +CD3X, 

Z s z , ^ • 

In subprogram AGIGAC «5 T'= P 
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aji'uiA 
b s h tl * 
‘H ns 
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' M .n s 

KaOMr 
Xs'O ( 
;isxi { 

As -A 
iU »\s 
Utils 

HO 

X’ M J 5 
X K ( ,) * 
X R n ) 

xp.n I 
>(Rt 3) 
^ F, ( ft ) 


• I . I .» It / \ t *i 

•J / R T ^ • .i 
Of.r-UT-TM 
h : h( b * 
r. ij F ( I' I 
f 1 1 

i ) 

A) 

?. J U . 1 

.-n 

~i I ( 1 ) 

s i » 0 

«' ■ I j • 1 ) *> y *Ul kl I #v I I 2 > 

0 > aiM g *n v.*V/( j ) »X 1 ( M 


aul. ( ‘•I 

= nl t •» > 
a » I ! 

a;: u ft ) 
u 


^ ♦•X'* r 
• A * X ft ( n 


stibP.nU’^ I :il; ArnsACi ST*Xn!iXRi 
ntfiR-sio'i xn3».xR(3) 

‘ib = ntui( ST) 

sbs.spn s"^) 

00 roi a, u » i 

shs-sb 
USX I ( u *st^ 

X P. t U =5; ’ ( n « c P - X I f 2 ) • 5b 

x R t «* ) 3 X M a 1 1 n f> ♦ \ 1 

X R t 3 I f- }J I ( \ ) 

MKTUrs 
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7 ilnttiiMDimtloYi nf nnnrditiaton and OoTiatitubnta of tho 
VVtvbbr ln_-t;ho Alaaoluto aYatom uf Ooordlftatoa 
Anoordlmt to the K'lnrami'iiR of Or'biHj 


1. FumrbimK accordin,‘T to tho Known olonionto of orKlt; 
a,e,t,«Mu^ aro dotorininod. 

f. dtructnro. Biibpro{':rani Kl/lBB 
Common units: /CHHl/j . 

3. Aoooss: OAII. BIABB (A, YA). 

4. Raw data: Aprav A.., containing valuoa of tbn orbital 

elements in the following oMer: a»ft. i,Ac^»U. All angles should 
be given in radians. 


5, Results; Array Y. , eontainlnxt values of X,Y,0 ,Yy ,V^ ,Yg 
in the absolute system of coordinates* 


b. AppliGvatinn of the area COMMON. Constants are used from 
the unit. aHtR/i (No, 5 table C.l). 

7. Algorithm: She following correlations are used: 

V - r (cos Q COSO. - sin Q sin a cos L ) , 

*/ =5 r (sin Q costt + cos Q sin u cost ) , 

Z ^ r sin u sin L , 
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iv = l^CcosJQcosu-sinQsin a cobi) u +sinQcost^ cost), 

Vy-\r Csin S^cosu4cosQsirmco8t) -V^(sinC5sir.u -cosQcoswcost), 
\i = Vf. sin u sin t 4 cos u sin t , 
where 

r=P/(l 4eco5 V) , t>att.-co, 

l(,=Cf4/P)'{^ e sinv, =(|4/P)''^V|4costo), P^a(hcf, 


I L .-product of the gravitational constant into the Earth's 
^ mass. 


8. lext. 

•liir'i'n.j’ ' ’IH ELAh'*<A»>$) 

) f ti I r I *, < A ) » •>< n ) , V < H ) 
nO'i'nN/C'.iH/r,’’ « 

f> s t 1 ' - I 1 , - A ( ? ) ** 6 < 2 > ) ' 

.10 1 

'(( J I = Ji! S ( A ( J *? M 
1 .1 ( - i. 1 - S ■ M ! A ( J ♦ '< 1 1 
,;(>)= M 1 o t ! ; 1 

; I o > c -r ( t, } *''.(• ( 5 ) 

V.) = S .-3 ’■ I r»H/ r- ) 

VH~" l“A < f > 

1 . *•;, { 2 ) *.) { T ) 

Vll-'' IV]; 


‘<s‘Vr. 

f's»M 1 '*w{ 1 » 

• M ,l < 

iH 3 »sj{ :») 

.n / » s-,M n ' 

DO / 

2 X < .J* ■» 1 rtw( J * '. ) •!»? S J J*5 ) •»/( I ) 
')0 ^ 

X < .» ♦■ ■? > = M ( .1 J • V P - V » .1 • 3 ) • *' 'j 

3 X ( .1 1 s X : .1) *» •< 

A ( 3 » 5A < 6) -A ( 1 » 
ftKTltR I 

i Eii:) 


5.8 D etermination of Orbital Elements and the /y29 

Position of Points on the Orbit According ■"” 

to Known Ooordinates and Constituents of 

tFe Velocity Vsclior in the Absolute System 

of Ooordinates 

(bOs-ab'sHTT 

1. Function: According to the known coordinates X,Y,Z and 

the constltutents of the velocity vector V„, VI*. V in the absolute 

X y z 

system of coordinates, the elements of orbit a, e, i, n, w, are 
determined as well as their position in the orbit u. 

ii. Structure. Subprogram ABSEL 
Common units: /CGR/^, /CFl/^. 

Access: CALD ABSEL (YA, A). 
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4. Raw data; 
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Array lAg, containing the values of X,Y,I3,V^,Vy ,7^. 

5. Results: Array Ag, which contains the values of orbital 
elements in the following order; (X,, ^ i ®he values of 
all angles are in radians. 


6. Use of the area COMMON . Constants are used from units 
/CGR/^,/0PI/5 (No. 5 of Table 2.1, No. 1 of Table 2.3). 

‘T Ai a=r/(2-k), ^ =rKVfJL , , 

7 , Algorithm. 9 a ' « 1/2 

r^^X.fY. + Z,, e » {1 -/t(2-*).cos*0.y , 

sin Q = (A'i^ 

Cm 

i = -- arct g , 0 < t ^ at , 



,c^=rV^-ZVy, c^=ZVy-XV^ , c^^XVy-YV^, 
c =(cj .cl+cf)*'* , 

kun Q cosQ 

9 


/ c,\ 

n =arctg t = orctg 

V ^ 2 / 


k cos*Q-l 


CD = u-'V 


ZC 

U *■ OfCLO y ^ . 

® YCrXCg, ’ 

[i. product of the gravitational constant into the Earth’s mass. / 30 
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8. Text. 

rUP'ROU'^ I 'lE AP-5 El<X»4) 

: 1 MENS I ON X { i ) » 4 U ) • C < '• ) I s ( 4 > 

' iMFj’SInN ",6) 

C0MM0N/c-'*>P/^4R . _ 

•CO».mON/C^ i/Pr*PP2,Pl2 

• -• « 

• w 

» »• «k 

> V S 0 

..ro ij-i,> 
p s X t j > " / Pi ) * ” 

. = X ( J P 

• i-v = x( .1 ) a'^a j 

f = SRRT.(P ) 

A ! r P. / '*♦ R 5 • *' 

_ s T R 5 
'r 12 J = 1P 
t ; 5 ? X t ,: ) / R 
7 ( I * M = y ! J + !V ) / V 

'■ f 1 ) = W f ? ) * - V P ) ® V P 5 

: ( 2 ) = V f ;' 5 «» V f 4 ) - V I ,\ ) • V f > 

‘ 2 f 3 ) = V P ) • " ( '■' > - V P ' » i, ' ^ i 

r = 0 ! '• - V- c ( 1 > ♦ c ( •? ) • c ( r ) • c P ) 4. r ( 

r r.= S'". ( '“Ml 

I’r ( C*‘'“ 1 P"'”' P A f I • , 1 

\y r ( •'. ) s c I > ) • X { 3 ) / ; *’ - r, ( r I • ’•' ( i ) / r ' > 

: i I’ ( 1 ’, = S R q ’• ( c ( 1 ) • r. 1 ) ♦ r ( 2 ) * r f 2 ) i 

3.9 Standard Array of Initial Cooditio-ns 
TB09-TRDAa)0) 
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' r 
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r, 
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1. Function. Subprogram TRDATO extracts from the initial 
conditions, which are given in the standard form information which 
is necessary for conducting navigational computations. 


Ihe initial conditions can be given in three forms: as the 

elements of orbit; as the coordinates and constituents of the 
velocity vector in the absolute (b) and G-reenwich (c) system of 
coordinates. We will designate as PN the array which contains 
the initial conditions in the standard form. There are 12 ele- 
ments in this array. In all three cases the first four elements 
and the second elements of that array contain the following 
parameters : 

PN(l)-the number of AES and the launching date in the form / 731 
of a decimal fraction, the first three numbers after the 
decimal point are the number of the AES, the following six 
numbers are the date (the year is indicated by the last two 
numbers) ; 

PN(2)-the current date in the form: C.EDMMGG, where DD is 

the day, MM the month, and GG the year; 
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PN(3)-tbe number of revolutions; 

PN(4)-the present time in the form O.HHMMSSDS, where HH is 
the hour, MM the minutes, and SSDS the seconds and fractions 
of seconds ; 

PN(ll)-is the ballistic coefficient (m^/kgsec^). 

Ihe rest of the elements of array PN are different for all 
three cases. 

In the case a):PN(5) is the Draconian period in minutes, 

PN(6)-the semiraajor axis of orbit (m). 

PN(7) -eccentricity , 

PN(8)-angle of inclination of the orbit from the equator (degrees), 
PH(9)-vertex of We orbit (degrees), 

PN( 10) -argument of perigee (degrees), 

PN( 12) -minimum altitude of the orbit (m) . 

In the case b) :PN(5)-PN(10) contain respectively values X,Y, 
PN(12)-2 

In the case c) :PN(5)-PN(10) contain values x,y,a,v„,v v„; 

PN(12)==1. X y a 

In both cases the coordinates are given in meters, and the 
constituents of the velocity vector in m/sec. 

In the v/ork results of subprogram TRDAIO the angles are trans- 
lated into radians, the time into seconds, and the calendar date 
into the RJD. The stellar time is determined for Greenwich mid- 
night on the current date, and the initial conditions are converted 
into all three systems of coordinates. /32 

Dimensional quantities are converted into the system of units 
given by the scale factors: EM, ESEC. 

Note: If EM=1, then the initial conditions remain in meters 

and seconds. 

2. Structure. Subprograms TRDAIO. 

Access to the peripheral subprograms: DAIDAY (BOl), HMSSEC(B03), 

Sa}TIME(BC5), AGIGA (B06) , BLABS (B07), ABSBI (B08). 

Common units :/CDEGR/jj^ , /CB3/^,/CC60/^. 

3. Access: 

GAIL IRDAIO (PN, NSP, DZ, RJD, NB, I, ID, SB, B, SO, A,YA, YG, EM, 

SSEC). 
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4. Raw data; Array RNnp, which contains the initial condi- 
tions in one of the three forma deacrihed in p, I 5 


SM, KSEC-scale factors. 

In the case of a, it is necessary to make the value of A(6)-u 
(most often the initial conditions are given in the point of the 
ascending vertex of orhit, where u»o)* 

5 . Results; NSP-the number of the satellite; 

D7.-the launching date in the form: O.RDMMGGr; 
RJD-the current date 'as the relative Julian date; 
HB-the number of revolutions; 

I-the current time in seconds; 

TD-the period in minutes; 

SB-the ballistic coefficient; 

B-minimum altitude of the orbit; 

SO-stellar time at Greenwich midnight on the 
current date; 

Array A,^-which contains the elements of orbit 

a.e, L\ n. 00, u 

( angl e a in rad ian s ) ; 

Array YA,— which contains X,Y,E,V ,V V„; 

Array YGg-which contains x,y ,ss ,vj,v' . 

6. Use of the area COMMON: Constants are used from the units: 

/GDBGR/^ ,/Ci53/^,/CCG0/^ (Nos. 2,5 and 8 of lable 2.3). 


7 . Text. 

.) .1 ■' h ' ' ; ■ I r. ” r A 7 c ( f II s n , e ?, I c ^ » '1^ ' ^ » s B » K • 5 e » “ I X A I X R * 
«• H ' •. > E 5 5 r, r. ) 

'>r>i*'S: I'll t » A I 6 ) . XA ( 6 ) 


an 

_ an a. in i/r nfiC!-' /oiinR 
anM.MC’i/cC'nu/r 

a » 1. 1 . .jAT'i A V ( mj t •» ) » D T ) 
a.;u. il,,SSliC ( nil) A ) ► T 1 
T = r / p s t- f. 

.TRKU'J ( i ) • f; A 
•J 5 PSTK 
n ?. S < T R - >! S P ) 

'•/-Bia.'iiaF.c 
SB = e>t( A .u 
;ih = i"i ( n' * . i 
r.Ai..i. S'' T 5 .(!;• I - ' » so ' 

= t X ..* 1 * • i 
! r ' ‘J - 1 ) w X 1 3 

X Til “ 4 “ - I 

I .) ) .1 * J, 5 /(• ". 

A ;; r. > j - n x = ■? 'i i j •*• ’ ) / v 

a A I 1. A i, 1 j 1 S n * T ( X t. I R I X A ) 

•'.OT,l 7 

^ ! r in - , 1 J ^ , R . 5 

: > n r. j .I ; , 1 


itAu ) =;'ii( j^a'V/em 
S XAt J*5 ) =PN ( 

CALI. A 0 I G A I S 0 • T . X A I 1 , X 'J 5 
CALI. AllSEL< Xi *A)- 
Tn=ii 
B = u 

GOTO n 

B no 

p A(.) ) = ?;;(, 

A ( X 1 r p,. ( 6 ' / E". 

A ( R 1 s Pi; ( 7 > 

A(p» = Ai6)/r)Enp 
To-p!,'< .•>) »af>D/FSF.a 
B = i'n( 

call EuASS<A,XA) 

Cali. Aaxr.Aisn.T.xA. 

K K T '•) n M 
r.M'» 
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3,10 Altitude of the ASS ovr»r the garth’s 

Surf aoe . the Geographic laii'bud.e 

and longitude, the Gedcentr^^^^ lati- 
tude of a Sub-Satellite 

TbIS^GBOGRG, HF.IG‘ril?. 6oi!l}ffn 

1. Function. Using the known Greenwich coordinates x,y,z, 
of the AI5S, the subprogram GSOGRO determines h— the altitude of 
the AES above the surface of the Earth's ellipsoid and,*<f » Jt — 
the geographical latitude and longitude of the sub-satellite 
point; subprogram HEIGHl determines only h; subprogram GOITIN 
determines i A the geocentric latitude and longitude of the 
sub-satellite point. For si.bprograms HEIGHT and GGITIN, x,y, and 
s can also be given in the absolute system of coordinates; in 
this case there is the right ascension of the AES. 

. Structure. Independent subprograms; 

GEOGRC, HEIGHT, GCITIH. 

Common units ; /CAE/g »/Q1*1/5* 

3. Access; 

OAlIi GECGRC (X,HG,AIT, A2iR, XRG), 

GAIL HEIGHT (X. HO), 

GAIL GGITIN (X', ALTO, AIN). 

4. Raw data; Array X3 , which contains the values of x, y, z ^ 
(for GEOGRG only in the Greenwich frame of axes). 

5. Results; 

HC-the altitude of the AES above the surface of the Earth's ellip- 
soid; AIT, AJN-the geographic altitude and longitude of the sub- 
satellite point; array XNG3 , which contains the cosines of the 
peripheral standard which is directed towards the Earth's ellipsoid; 
AlTG-the geocentric latitude of the sub-satellite point. 

6. Use of the area COMEON. 

In subprogram HEIGHT constants are used from the units 
GCMMON/CAE/2 ,/0ASI.A,,/GCI^V'i/GP 1/^ (Nos. 11, 12 and 13 of Table 2.1 
and No. 1 of Table 2.3), in '^subprogram HEIGHT constants are used 
from units; /GAE/_ ,/CAE2 /t and subprogram GGITIN from unit 
/0P1/3. 

7. Algorithm 

h = r - ttc t a^ctz^/r^ , 

<p = arctg [z/r^ ( I -a/} , -ji /2 < <f>< Ji/2 , 

X =Arctp ( y/x)y 0 < A-c 251 , 


where f + 2 “ , 


) 


2\i/2 
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ag,a -is the semiraajor axis and the contraction of the Earth' 
whole ellipsoid. 


cPoc*arctg(«;/r,) -the geocentric latitude. 

where ^ ^ the cosines of the peripheral standards directing 

towards the Earth's ellipsoid, 

8. Texts. 


St«»RoUTlNE GEt)T;RCiX»HC* aLTi ALN» xe ) 

COMMON/CCL2/CtZ 

comhon/cae/ae# Al 

COMMON/CP I/pj »Pl02»Pl2. 

common /CAEl/AEl 
dimension UC2) .cU(2) »X<6» 
Rl2sxU>*XU)*Xl^ )*X< 


il*X< 3)*X(3) 

R2=ftl2+W 

HCsSaRxCRZI-ftE-’-AEt^W/RZ 

CUil)=j;.<i» 
CU<2j=CLZ*SQRTCRa2> 
l’|sS(lRTOV*CU{ 2 )*cU( 2) ) 

DO A jsl,2 
ir<tu(j))2*3#2 
3 U( J)rsiGN<Pl02 * X< ) 

• goto a 

2 U( J’> = A7AM‘XJj^1)/CUCj)) 
A :<E‘ J»*CLZ/W 

1 continue 

XF.‘ A)=x<3)/W 
ALT*0( 2) 
alnsih i> 
lr<CU( 1) )S»6,6 

3 AlNsALH-. p! 

6 Ift AI.N) 7,b,a 

7 AlN = ALm'*PI2 

fl return ■ 

End 


SunKnUT INE he I CMT< V » HC ) 
-dimension y(A> 

COMMOli /CAE/AF,»AL 
COMMON/CAKL/AF.!- 
MSM( 3)#yt 3 } 

UO = .SRRT ( R2 ) -Ac'*’AeL*W/R2 ‘ 
Return 

END 


suRP.nUT I nf nci tln ‘ V » alt » aln ) 

COMMON /CP I /p T , 9 I [)2 , P I 2 
DiMK'iS: ON N(3) 

RL--Sr,RT<y< :)-^V(2J«y<2) ) 

A L T = P 1 D 2 
1 F « R 1) : • 1 • 2 

2 eLT = A7/.»:(V<3)/Rl) 

1 ai,n = at,\?:;m y t 2 ) , y ( 1 ) ) 

IflALNliifc./V 

3 AlN = ALI:*P12 

4 HETilRii 
END 
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OHAPfER 4 

POSITION CP IHE MOON AND EHB SUN 
(INDEX 0) 


/36 


4.1 Position of the Moot) 

(OoI-sIM a)' 

1. Punctlon: The subprogram determines the poijition of the 

moon in the absolute geographical system of coordinates for the 
moment of time which is given by the relative Julian date and 
Moscow time in seconds. 


2. Structure. Subprogram SEIENA. , , , , , , 

Oommon units: /CRS/3_,/OSDAy/;^,/OT3/3^,/CDSJS/3^,/OPl/3. 

3. Access: CADI SEIENA (RJD, T, XS, RS). 


4. Raw data: RJD-relative Julian date; 

T- Moscow time in seconds. 


5. Results: Array XS^, which contains the direction cosines 

of the radius vector of the'^Moon in the absolute system of coordi- 
nates; RS — the module of the Moon’s radius vector. 

6. Use of the area COMMON. Constants are used from the 
units : 

/ORB/i,/CSDAY/i,/CT3/i,/CDSJS/i,/CP 1/3 (No. 9 of Table 
2.1, Nos. 1,2,4, of Table 2.2, No. 1 of Table 2.3)* 

7. Algorithm; 


X° = COSOfij C0S &J= COSPg COSA.g^ J 

y® = sin cos&4= cos pg s‘m Tt^cos e - sin sine , 

zl - sin =cos Pepsin Ag sine + sin fg cos e . 


In 

in [ 71 , which 
an accuracy to 


accordance with Brown's theory of approximation, pregefited 
guarantees a prediction of the Moon's position with 
30", the following correlations are used; 


G » 2 3 ®27 'Oa^26 - 4 6 ", 8 4 9 r- O", 00 5 9 r O' 0.0 1$ I T 
^ - 296 ° 06' 1 6" 99 + 4 7 7 f 9 S ® 50 ' 5 6 '', 7 9 r + 3 3 "09 r O", 0 5 < 8 r ^ j 
28' 33'; 04 + 69999 ^ 02 ' 99 ", iO T-Oj 94 T o’o^20T^) 
H®l9'03",20+483202®0i'30'547-ir,56r*- 0*0012 7^®; 
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+^-^ 5267 ^ 06 ' 5l^<87'-5'j!7T^+0.0068T^; 07 

Ag = 270‘’26'02',99 +-48 1 267° 52' ^%HT-A’0&T^->-O.O.ObST^->- 
+ 22639* 58 sine --!i586’,/<38 sin (8 -2D) +2369' 899 sin2D + 

+ 769”, 021 S)n2£ - 668", 94^1 sine'-^ii",6l4 sin2f- 
-2H',658 sin (2e-2D) -206* 219 sin (8 +e'-2D)+i9i,95^«n(e+2D) 

-i65*,35l sin(e-2D) +i47* 878 sin (e-e')-i2<i,785 sinD- 
-I09*,804 sin(8 +8')-55*, 17-A »in(2#’-2D) -49"l00sin(8+2/'’) + 

* t 

+ 39",5&2 sin(e-2e’)-38”^28sin(e-4D)+36",i2/i sin 38 
-30'773 sin (28-4D)-28* 5li sin (8-8'-2D)-24',45i sin(8 '+2D); 

=iS46i*480 sitif +iOiO"l8sin(8+/!') -999'j695sin('A’-8)- 
-623'.65S sin(8’-2D) +l99",485sin(F+2D-8)-i66",577sin(8+8’-2D) 

+ H7*262sin(r + 2D) +6i ,9i3 sin(28+8’)“33*359®^^(^“2D“8)'- 

:n:, =3422",7 +iS6'',539£ c.os 8 + 34”, 3H? cos(8-2D)+28,2333cos2i) 

O f 

+ I0'i657O-os 28 +3'086lcos(8+2D) + i',9202c.os88'-2D) + 

+ i",4455cos( 8 +e'- 2D) + l'i542cos(8-8')-o'9752cosD - 
-0'9502cos(-8+8')-0"7i36cosf8-2F)+0"62i9Cos'3e + 

+ 0"6008 cos(8-4D); 

Z„ = 206264",81«3/:Hg, 

where p is the equatorial radius of he Earth. 

3 

T-RJD/Wsb , 

RJD relative Juliao date in days, 

T-the same interval of time in Julian centuries 
■C -the median anomaly of the Moon, 

1 -the median anomaly of the Sun, 
f -the median argument of the Moon’s latitude, 
p -the difference between the median longitudes of the Moon 
, and Sun, 

■ -the Moon's longitude, 
iPe -the Moon's declination, 

^ -the Moon's parallax 

Jig 
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8. Text. 

SUttRoUriME SElEM*<OT»T»SiR) 

OIHEmSION SO) 17) 

COMMOII /C7S/T3 
COMMON/CPf/PI iPI92*Pt2 
COMMON/COSJS/PSJS 
COHMON/C5DAV/50AV 

commom/crg/rc 

TRS0T/D5J’S*<T-T3)^SD|'^/0SJ5 

5U»sM. 6773127636-0. TP -♦2860 40 0 7 2E-7>*TR 

- -.227110969e->)*tR*, 6C931975S3 
W»6325 ,69l*TR 

W = W-*J*P!2 

A(4)s{<,2»U33487E-6.TP*.l6ti42fc8*7p-3)»TB 
♦ • •,i03»697lE-3)*'TR*-6.1*6000347*W 

US626. 30l*TR 
NsW/pta 

A{27)sn-.9817764.17E-7»YR-,261799 38 8C-S)*TR 

+ r, 04> a 7c20iE"3) •TP. *6.25656 3 776-1) 

i|s»t33 , C66* TR 
uaii/p :r. 

U£i>!-.)«r ! 2 

56:77641 7 56 0 > F. “4J*7R 

- * , ?.9 5a9’*.‘i3’»E-3 ).TR>-, 1963 65: 540-Vi 

us 7 ^71 , 37 7 ‘TP. 
fisrt/ p I 
wsw-Kis-pTn 

A(2>s<<,32y6733e3E-7*TR-,?S06^9673E-*)«TP 
• • , 1967u89207E*3).TR*6 . 121523V44i*W ‘ 

4{3)=2,.a<2) 

7,<?)=.‘»,«A<4) 

00 i J“4,7 

1 A ( J - 2 > s 4 { J V - A f3 ) 

‘00 2 J-; 4,6 .2 

4(J*")=fl( J>*A<27) 

2 A{4*6>3A{.i)-A(27) 

0 J 3 J 3 3 . > 

l.s2»( j-l) »14 

A<L»saU7)*A(j) 

3 A(U‘l)=4{l7)-A< J) 

A(ll»3Al27)-A(3> 

AU>s3*AM*> 

A(24}s2.a( 17) 

A(25)sa(24>-a (3) 

A(26)=A(4)*A<29) 

A ( 1 Q > = A < 4 » - A ( 2 4 ) 

A {26 ) =A( 15 ) (17) 

A(29>sa(6)‘A{ 17) 

A(16)3,-\(4)*A(3) 

A { 30 ) =A ( 16 I -A ( 4) 

AC3n=A(27)«-A (3) 



/S { .'2 > *A t 21 ) - A f 3'^ 

00* 

A o«.i*x)snosuu>' 

0C>*. 

5 A { J> s s * N r A < J > 1 l 

0 ( ^ . 2 1 2 ? 6 0 6K - (3 ( A ) • , 5 0 n I i7 0 :5g* S*-8 { 2 5 

-•-»3A37^30‘‘3E**S*6^ H > «*i. A00ft»9?»6|j«'9-’I5 { A?. J -• A F 'ST'»{J ?a* e *» 1 

*. **3 59 57iO>Xn-5«8fJl3)'‘ . 7 <! 0 7 9B ‘*3^5 i X.l- « » ,i. * 

^,l40 6lf3?g*4«fi( 17U.A920A7riAA|j-4.ni /.I 

* * f J A6347B X ( 7 I .V .9()^3T0A7j.B-3-ft { » J Jft»9 "TX7P«»« j. 

Wscsi99,7oy^7R 

tiT.l^/9:2 

W s W •» *1 • n I 2 

rtf ! • s". L'/S«;->‘?3E-3^fi(3X)-, l3ft?2»2?9F-3*/» UA5 

" - . i * - * Af t i 9 } * , 1 7513 AftO* A ! J ) - , XH6?fS»- J4BXH-3-A I B I 

* •>..‘■l63^.'j«*4^c-3•At^3)••.2if65097E-3#-A(^e)-.?.67A9lXf!-3»*A^^!5) 

, « , ‘ i 4 j;- 3^A< XI fsO' 9 7A75jg.3^4 (Pj 7.t ^ 93 2?7A ?- 3* A < 12 ? 

- - . ' ■ ’ ’r.--^^. ,;' I5^^.«^f35i9?.53c-3«A( lc,J..0 9«^F"F9f‘5F.«j#A([ 13) 

- -, • - ‘.*Uo9<.r-2«t' ’)-.;.99S56?99B-2*A<2A)-,3;;A3lAi{93E*.i»*A(27) 

^ *.. s- • V v J2f. , p,. ! 5)* , J»A ( 3) - . 222356 7e9g.jt*A( 6> 

4 ' 1*A < 4 ) *', { , 32e673303E-7#TR- , 197*P * 932I:-A ) • tb 

* * » ,'. 4 fj 'p fc ^ V 7 !■ - i } «. - i , 7 j 9 9 ^ 6 5 7 Q g *11 
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f 3 v = - , li.5V3«s33tF - 3».*. ( 2B ) - . 1539 9l37E~3vfi( ?3> 

- -.l‘>i''23s»vAf;-JM<J2}^. 30 016269AE*3*A(22)**36e9U22l9E-3*A(lfl) 

- -♦'‘*’5'-'Af'»*$r-3'*A<29)*,967no572r:-3»fi(30)-iSc233793lE-2«A(l9) 

- -*‘-fttaSA‘-'ii3c-2.6(2i)i.,4697490£/. t*-r»4f2n5*.8 03337eCfBP-J*A(17) 
’) r 5 . n ', » 9 » 2 

h ! J • •. } = ^ '- <> f n t p ) ) 

t h ' vM s< ; > ( H ( ) ) 

A 9 » £ B 1 6 ) ^ > 

. J. ( J • £ : ' - ' • rt ( n * f. { 5 ) • P ( 2 ) 

h / 2 ' s f’t ; ) * *• 5 \ " f* ( '5 ) • fe ! 1 ) 

S ( I ) r .- ( *< ) • ( 5 » 

k 5 H E / • £ 1 

H gum •, 

KtlO 


4. 2 Position of the Sun 
(C02-SUN) 

1. Function. The subprogram determines the position of the sun 
in the absolute geographic frame of axes for the moment of time, which 
is given as the RJD and as Moscow time in seconds. 

2. Structure. Subprogram SUN. 

Common units : /CAED/^ ycSDAY/^ ^/CT3/^ ,/CDSJS/^ ,/CPl/ 3 ,/BIECL /2 
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3. Access; CALL SUN (R«3D, T, RS, AS, BS, XS). 

» 

$, Raw data: RJD-relative Julian date; 

T-Moscow time in seconds 

5, Results: RS-the module of the Sun's radius vector; AS, BS 

the corresponding right ascension and declination of the Sun; array 
XSj-di recting cosines of the Sun's radius-vector. 

6. Use of the area COMMON. Constants are used from the units 
/CAED/^,/CSDAY/^,/CT3/^yCDSJS/^,/CPl/3 (No. 1 ofTable 2.1, Nos. 1 
4, of fable 2.2, No. 1 of fable 2.3). Into the unit COMMON/BIECL/g 
dispatched the values co5£,5ine, which are defined below. 

= 2 > 9 ® 10 * 59 V? - *9^4 ®0S' 3 1", 2 7 0 *, 006 T ^ ; 

£ = 2,5®27'08'26-46>457’-o’0059T^+0,00l6!r^+9’',21co5S^; 
e^^0,0\6y^lOA -0, 00004 isr -0,000000! 26f ", 


7. Algorithm 

= ArclgCsinXQCOse/cosA.Q) (0,06! 164 1 ^Ay)'-20 , 496 ; 
60=arcig(sinA.^5ine/(co^A.0-t sin^X^cos\y^^)-20' 496 sins coaa^', 
A.Q = X -»-2eQ3in (K~^)+Va sin 2(X -Jt ) ; 

'"o' “oil + 2M^)/i]\ 

it -28l'’l3’l9'+6189’,03r+l',63r*+0'l2T^; 

Aw = -17,23 sinCJ ; «0*l,OOOOOO23a,c. ; la.e. = 1A9600000 km 

tA =A-3t .^=279‘’4r48'04 + l29602768',13r + l',0697" 

where T=RJD/36525 

RJD-relative Julian date in days, 

T-the same interval of time in Julian centuries, 

X© scosa© cos 6© , 

r© = si not© cos 6© , 

Z© *sin 




2, 

are 
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CHAHOSR 5 

MODELS OP LHE’eARIDH'S ATMOSPHERE 
(INDEX D) 


M2 


5.1 The Standard Plve-Iayer Model of the 
Earth *s Atmosphere ’ 

(SOllRg) 

1. Punction: Subprogram RO calculators values for the 

density of the Earth's atmosphere p at the point given by the 
altitude above the Earth's surface. 

2. Structure. Subprogram RO. 

Ooramon units: /GHA/gQ. 

5. Access: RO (HC^P) . 

4. Raw data: HC-altitude above the Earth's surface. 

5. Results: P-the atmospheric density. 

6. Use of the area COMMON. Constants are used from the 
unit /CHA/go (Table 2.4). 

7. Algorithm 

where h is the altitude over the Earth's surface, i ' ”2t 

are the coefficients depending on the altitude h (Table 5«1). 

TABLE 5.1 Values of the Coefficients 


\ » 
1 ^ 

*>• 

/A 

A 

Ai(Kr ■ 


k2i,(Ui'') 

Ju 

t 

IOOCOO< A <150000 

0,4141 

0,1469 10“^ 

0,1767 I0"“ 


■ 150000 «/?< 300000 

0,2173 10*^ 

0,8004 

0,3734 lO"^ 

! 3 

3000000<;/l< 600000 

0,4861 10”^^ 

0,7111 lO’^* 

0,1547 lO"^ 

1 4 

60000Q0^1< 900000 

0,B904 

0,1831 10"*^^ 

0,9275 10“^ 

1 5 

900000^/t 

0,649? 10"^^ 

0 

0,S540 10*® 



ORIGINAL PAGE IS 
OF POOR QUALITY 


8. Text. 


■> vl K 0 J H P, t '< ** 1 5^ ) 

C D/I.M'’ 's / r >1 .‘t / H ^ i ; 5 ) , (j ^ { j) ) ^ j 5 j 

Dn*i Js2 . 5 

I F MC”::.v( j > M , 1 , 1 

1 I = I*l 

2 CON T I -lUK 

3 ( I ) 

p = aa( : )*EXP( (ni ' r »*w-Q2{ i > >#v»i 
HF. TIIH U 


5,2. Model of the Earth’s Atmosphere, with a 

the Influence of Solar 
ITadio Emission at Wavelength 10.7 cm 
of Geoma, Kinetic Disturhanoe, Daily and 
Semiantiual Effects 

'mr-vwA7'''mwsi^ ^^^^ 

1. EuTictioT]. According to the given positions of the point 
and Sun, and according to the values of the indexes Ei|o-<7 (the 
intensity of solar radio emission at wave 10,7 cm) and 'ap 
(three hour indexes of the geomagnetic disturbance), the subprogram 
PENS determines the atmospheric density at the given moment of time. 

Coefficients of the model of the atmosphere are chosen with 
the help of subprogram VKMA depending on P , the mean level of 
solar activity. 


2. Structure. The packet of subprograms* 


VKMA, DENS. 

Peripheral devices; printer 
Common units; /AKOEP/^q,/SYEAR/^q. 

3. Access to VKMA; CALL VKMA (F). 

Raw data; P-the mean value of solar activity. 

Possible values of P; 75, 100, 125i 150. If P differs from the 
indicated values, then a stop should occur, whereupon the printer 
reads; "P incorrectly given." 

Results; Into the unit COMMON/OKOSP/A(50) , are sent values 
of coefficients of the model accordine: to Table 5.2, and into the 
unit COMMO'N/SYEAR/I^.§r is the set of numerical corrections on the 
r -semiannual effect accordinp; to Table 5.3. 


^Author of the subprogram M.I. Voyskovskiy. 
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4. Access to DSNS: 

fiAlL DENS (HO,X,Y,Z,SUN,AP,n 07 ,D,I,RO). 

Raw data; HC-altitude above the Earth's surface in km; ^/44 
X,Y,Z, -directing cosines of the radius-vectnr in the 
absolute system of coordinates; the array SUN, which 
contains the values a 0 ,o© (right ascension and declina- 
tion of the sun in radians); 

AP-th® value a p at the moment of time where 

t is Moscow time (for isolating the calculation of the 
influence of geomagnetic disturbance it is sufficient 
to let AP< 0.5); 


F 10 ,^-the P^Q , at the moment of time t-AXp where 

t is Moscov/ ?ime (for isolating the calculation of the 
influence of solar activity it is sufficient to let ^10 

< 0.5); 


7 


D-the date and Moscow time in the form of the number 
of days from the beginning of the year (for isolating the 
calculation of the semiannual effect it is sufficient 
to let D < O) ; 


I-the parameter controlling the calculation of the 
daily effect (DE); 

if I <0, then DR is considered memberless with the 
coefficient Og; 

if 1=0, then DS is not considered; 

if I> 0, then DE is considered in its entirety. 


Results; Array ROg, which contains values of five factors, 
each of which takes into consideration its. effect in the modes 
of atmospheric density (see p. 6 ) and the value of the density: 








5. Use of the area COMMON; 

In subprogram DENS coefficients of the atmosphere model and , 
corrections for semiannual effects, whose values are assigned by 
a preliminary accessing to subprogram VKMA, are used from the com- 
mon units /QKOBE/^q and /SYBAP./^g. 

6 . Algorithm 

I'he density p is calculated according to the formula: 
p - py' k ^ kry , where fi_, is the nightly vertical 
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c.'oss section of the atmospheric density. 

k, -is the factor which allows for the influence of measure- 
ments of the intensity of solar radio emission S' at wavelength 
10.7 cm relative to certain mean levels of solar radio emission. 

,k2-the factor which allows for daily effects in the disper- 
sion of the atmosphere, 

k^-correction for the semiannual effect. 

kij-the factor which allows for the correlation between 
changes in the atmospheric density and geomagnetic disturbances. 

'V " ' +(C; - Cj, h +o_ exp (-(h +c^)V4|(cos”’' ij j/2 

^a= 1 + (Aj +Aj/;),4 (d ) ; i In {ap/d, j 

where cG& shti 8 q + cos6q(X‘'<ios +y^s-m vj > 

co£, =-Z*’sin 6g, + cos6o(X‘'cos Yj +V‘siri Vj)., 

Y,=cx^+<p^ , V2=ot^-i-<P2 , 

// -the altitude above the Earth’s surface. 

, t .o . c-the directing cosines of the radius-vector of a point 
^ absolute system of coordinates, 

^■'o < -the^right ascension and declination of the Sun. 

d -the date and time in the form of the number of days, counted 
from the beginning of the year, 

i (cf ) -correction for the semiannual effect (the volumes of are 

shown in Table 5.3, with 10 day increments of time, the inter- 
mediate moments of time A(d) exists as a linear interpolation), 

f -the intensity of solar radio emission at 10.7 cm, /46 

Ctfi -the three hour index of geomagnetic disturbance. 

The values E and^/> should be at the moment of time, 
respectively: ^ ~t-Ax wheie t 

' ’j Moscow time, and ^ A^ap '^^® "time lags" of 

changes in the density of the atmosphere which correlate to a 
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correspoTidiTig change of the quantities F andU2p . Ihe values 

At and AT^are shown in lahle 5.2, and also are A^f* * , 

and,, At ap , the duration of intervals in the course of which the 

quantities F and ap are maintained "by the constants and talDular 
values which are equal to them. 


In table 5.2 are also the coefficients of the model: CtijCtjiCZg,- 
which depend on the mean level of solar radio 


i>, '^ 2 » ' ■ • » 
emission F 


ap\ which depend on the 

0 * 


lable 
on 


le 5.2 Values of coefficients of the atmospheric models depending 
the mean level of solar radio emission F^IO'^'^ W.m"‘^Hz"’^(\;nits 
BF/Q50)). 


i r 

Desi g- 
lati ons 

Q(i) j 

( 2(0 


QCO 

I 

¥7~ 

75 

100 

125 

150 

Z 

CL. 

-14,030 

-15,095 

-17,028 

-16,072 

3 

a 

0,9108 

0,8229 

0,7198 

0,7155 

4 

c 

a. 

59,77 

68,92 

93,36 

70,33 

5 

o 

f 

~0 ,630 

-0,750 

-0,710 

-0,765 

6 

62 

0,00506 

0,00560 

0,00562 

0,00571 

7 

■ ^ 1 

0,130 

0,172 

-0,274 

-0,247 

. 8 

el 

0,00014 

0,00217 

0,00257 

0,00199 

9 


3.733 

3,784 

4,048 

4,698 

10 


-507 ,95 

-566, il 

-632,63 

-7v77,58 

• II 


189,85 

200,97 

230,76 

278,35 

12 

^6 

-0,041 

-0,047 

—0 ,038 

- 0,012 

13 


4.2 

4.1 

4,4 

5,2 

14 

tn,) 

6,0 

6,0 

5,9 

5.9 

15 


37“, 4 

34 °.'' 

34°, 5 

33 °, 8 

16 

■ i 

<'-P2 

325°, 9 

0 

0 

00 

n 

CO 

308 ° ,0 

322°,2 

17 

a; 

1-0,602 

-0,526 

-0,513 

-0,607 

■ 18 

Ai 

0,00669 

0,00636 

0,00631 

0,00670 

19 


"■0 , a.32 

’“-O.B'- 

- 0 , 12 a 

‘-O.IIb 

20 

\ 

GO 

0 

— 1 
0 
0 

0 

0,00x04 0,00095 

0,00089 

21 

22 

23 

24 

A v^' 

Up 

^^'CLp 

39,^ 

10^,5 

4 " 

ih 

ooh 

ic/Js 

4^ 

39 " 

lo'js 

V'' 

• 39 ^ 
icfjr 

I" 

bi 

- P 
Ctry 
S 

2 

2 

3 

4 
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lable 5.3 


Corrections for the sewlaoo'ual effect in dependence 
on, d and the number of days since the beginning of 
the year (unit G0rW0N/SI-5AR/P(58)) . * 


K 

_Y* 

i 

/>!•/) '■■■A Id) 


i 

d 


i ^ 

i 0 

i 

1 -0,CoV 


20 

190 

-0,172 

0 

10 

-0,033 


21 

200 

-0,180 

3 

/ ■ 

f-*0 , 094 


22 

210 

-0,183 

4 

30 

-0,088 


23 

220 • 

-0,179 

5 

40 

—0 , Oir/3 


24 

230 

-0,163 

6 

50 

-0,0'05 


{ C'5 

240 

-0,133 

7 

50 

0,039 


2S 

250 

-0,086 

! 8 
> 

70 

0,090 


Zi' 

260 

-0,018 

. Q 

I ^ 

O'J 

C,I23 


28 

270 

0,059 

10 

SC 

0*133 


23 

280 

0,123 

XI 

100 

0,126 


30 

290 

0,161 

n 

lie 

0,099 


SI 

300 

0,170 

13 

120 

0,059 


32 

310 

0,156 

14 

130 

0,017 


33 

320 

* 0,119 

15 

140 

-0*027 


34 

330 

0,073 

16 ^ 

'150 

-0,065 


35 

340 

0,027 ' 

r? i 

160 ; 

-0,103 


36 

350 

-0,023- 

IS 

170 1 

' -0,136 


37 

360 

-0,055 

Trv 

XD * 

! 

1 CD 

I O 

4 

« 

1 

1 

-0,158 
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Yakkia-72 Model of the Sarth'a Atmosphere 

rm-ismTimKK: -(^mTirwui.w — — 


/49 


1. I’upctlop, Acordipg to the givep locatloo of a point, o.c 
the Spp, the values of % 07 » ^ given moment in time 

determine the temperature 'at a given point, the exospheric temp- 
erature, the densiT^y of the atmosphere, the number of molecules 
of nitrogen, oxygen, argon, helium and hydrogen in a single voume, 
and also the mean molecular weight of constituents of the atmosphere 
(presuming that it consis’cs only of the elements listed above). 


2. Structure. Subprogram ADEN. 
Internal access; AMBAR, GRAY, TIOCAI. 


5. Access: 

CAIL ADEN (AMJD, SUN, SAT, GEO, TEMP, ADION, AMHW, RHO). 

4 . Raw data ; 

AMJD-date & time in modified Julian days and fractions 

of days; 

SUN-the array which consists of two elements: 

SUN^lj-the direct ascension of the Sun in radians, 

SUN(2)-the Sun's declination in radians; 

SAT-the array consisting of three elements: 

SAT(l)“tbe longitude of a point in radians. 

SAT(2 )-geocentric latitude of a point in radians; 

SAT(3)-the altitude of the point in kilometers; 

GEO-an array consisting of three elements: 

GE0(l)-the value of the index time lag-1.7 davs), 

GEO(2)~the value of the geom^etffi^^ index A , 
considering ^ /50 

that 3 = 2 . 667 , 3 = 3.333 and so forth, the time delay 

= 0.279 days. 

5 . Results: TEMP-an array consisting of two elements TEMP 

(1)— the exospheric temperature over a given point (in Kelvin 
degrees); TEMP (2)— the temperature at a given point (in Kelvin 
degrees) . 


ALION-an array consisting of six elements: 

ALIONt 1) -common logarithm of the number of nitrogen 
molecules in w , 

AII0N(2)-coramon logarithm of the number of oxygen 
molecules in ; 

ALI0N( 3 ) -common logarithm of the number of oxygen atoms 
in M^; 

^The subprograms used are from [ 6 ] and were tested and modified 
for a high speed electronic computer BESM-6 by Ye. Ye. Ryasanova. 



i 


AiION( 4) -common logarithm of the number of argon molecules 
in ^ 

ALION(5)-coramo^ logarithm of the number of helium 

-common logarithm of the number of hydrogen mole- 
miles in M?i 


cules in M- r 
,10 . . ^ . 

RHO-density (in Kg/m^). 


AMHV/-mean molecular ^eight; 


6. Due to a lack of space there can be no description of 
the Yakkla-72 atmospheric model or of the texts of subprograms 
which were published in f63. 
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0HAM3R 6 

ANOKAIIES Of THE EARTH'S GRAVITATIONAL- FIELD 

(INDEX 




151 


6 . 1 The Accglsratlott Vactor. DRtermlTied 
Sy the Inf luen c e of Angmai i e b of 
Earth *s Gravitational Field 
TFOI-MGa' 7 5 EgF , DsGIp : 

1. Futiction. For a point in space fixed Py the Greenwich 
coordinates; x,y,si, one determines the constituents -Ao atid 
of the acceleration vector which is specified ^ tne®5Hflu- 

ence o? anomalies of the Barth's gravitational fi«ld 1.2 : 

— the radial constituent of vector .Ag (project 
on the radius-vector of a point with a minus sign^, .i-C 
the merdional constituent (the meridian directed to the north), 

^8c -the projection of on the perpendicular to the 

plane of the meridian (directed to the east). 


Suphrogram DEG2 determines the acceleration vector which is 
determined by the influence only of zonal harmonics of anomalies 
of the Barth's gravitational field. Subprogram DEG5 determines 
the acceleration vector, allowing for the influence of zonal, 
tesseral and sectorial harmonics of the Earth's gravitational 
field. Subprogram D3G4 determines the acceloration vector allowing 
for the influence of only the harmonics 22, 50 and 40. 

2. Structure. The packet of subprograms; 

DEG2 , DEG5, DEG4. 

Common units; /RAD/g ./CRZ/.j^,/BcnNGR/5^g,/0A2B/4. 

5. Access; CALL DEG2 (X,DG, NM), 

CALI D3G3 (XG,DG,NM), 

CALL DEG 4 (XG,DG). 

4. Raw data. Array XG^, which contains x,y,z — the Greenwich 
coordinates of a point? ^ 

NM-the number of harmonics considered (NMl 22). 


Note. When accessing subprogram DEG, one may use coordinates of /52 
a point both in the Greenwich and absolute system of coordinates. 

5. Results. Array DG, which contains the components: 

of the acceleration vector. 

6. Use of the area COMMON: 

Before accessing any of the subprograms it is necessary to 

put the values + y^‘ -f z^)^' and Rl=(x^+y^)'"' into the unit 

OOMMON/RAD/R,R1. 




1 

Author — Ye. Ye. Ryazanova. 
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In the unit a0MM0N/CR2/Rg thi valu© R^ should be fixed. 

(No. 10 of Sable 2.1). 

In order to aecure the work of subprograms D302 and DE03, the 
•values of coeffioients of the resolution ratio of the Barth's 
gravitational field (according tn Sable 2.5) by means of a pre- 
liminary accessing of subprogram C0NGR(A02). In subprogram DBS4- 
an additional variant of coefficient values of the Barth's gravi-^ 
tational field is used from unit C0MMCN/0A22 j/^ (No. 4 of lable 2.1). 


7. Algorithm, 

a) calculation of zonal harmonics 


b) calculation of zonal, tesseral, and sectoral harmonics; 
^ntcef rt*i ^ 

Ap..=l/r 2 (''?/>*) 

2(««mC05/»t+p„„ sin/nt)(a,mj(sint)-«tg\|)/^n4 

ml 07=0 

m2 /w»d 


x,y,z-the Greenwich coordinates of a point, 

Ap, t-the geocentric latitude ana longitude of a point, 

L -y/r^ , cos L -x/r^ , sin =K/r , cos-^ 

-Legendre's joined functions 

/. , ,,when n :> m. 
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P„„ =(2”’-‘>eoS'l’^n-l,/)-l > 

an n m, Pnm-0 

sin m I ^ sin -1)1 cost ^cos(m-i)L &iii t ^ 
cos m L = -sm(/T)-t )l Sin L +cos(m-i iL t , 

P^o 1 ^(0 « n /^ ; shiS|)~{)/2, 

Pj , = 3 sill's!* cos ^ 3 ^ r ® j 3 cos* \s , 


il8 


I 
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-anomaly coefficients ^ % 

- -the greatest number of ^AicmniQB cotieidtred* 


C) calculation of only the first four harfonl 

xg A Sm 


cs 


^e,-4(fK7- 


•Pr^**** 


«2 


+18 Pi! ■'■ 


r(‘ 

o f tL 




<f+ 


3 2' ,V., 

+ " ( -"-g ~ *t 


/£Y 

■ly) i 6 V-'- 
.2 1* » • 

r\» . 4- 9 ^ * 

2 / ' V® 

<!>U ^ 2 


•" S / / 


^ 2 7\ P 2 r 

*o -l/£f il(-!2 ^^>121+6-^- f>®~ 

^§/:“ r\r) r \ ' 


/ J-- / 


V' 
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8. Texts. 
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OP'lENSIf/H K(31 »CS<5>,PWCi5i V i AUiri 

COMMOK/fjCOHeR/Ar^Mf 2T55 » 3Fv*5'l£^^?# 

COMHON/RaO/RiiR 

COMHON/CRZ/Rb 

aU»eA|jM(<») 

AUtsAUM(S) 

AOtcANMI I 
fi(A»sAUMUV) • 

A{6)fA«M(5i> • 

A< y)=A«M{ <^55 
A(tJ )sflf{«< 

A< VlCflfJMl 6.4 ‘ 

,AU8)«AKM< ("w 5 
A(XI> AANM« ’ 

Atl2}AAMM< iOSj * 

.A< X5)sanm< liBJ * 

AU4) ?ANMU54i • 

• » WH O • . » »'•» • ' • 

16 ) =Aii:i< 16':' > 

.'. f \ 7)=AN'.Mt l-“. 'i) 
i ( IB) 200) 

A { 49 ) =ArJM<,22'» ). 

r,l~Rt;/P. 

C2=X( 3 ) /n 
p f|U ( 1 > =1,2 

:ih|U(2) = (3,U*c?.*c2-l,Q)/2,0 
? H I U ) = R I . 

PMl i > spill ‘ I ) /R • / 

3'f’- ' 2 ) S.J ..i)*c2*PfU <1 ) / 

.)G^l) = w.O / 

) G I ?. ) = 0 . 0 

11 = 0 , 0 

DO lO . !I = 3 • f'M- 
U s N 

0 3=(n,0'''.J"1.3)*C?. 

Pr|li « -5 ) = < C3*P*;0 ( 0 ) - I U - 1 . 0 ) ‘PIIO < I ) ) /I’ 

P:jil3) = (,:3*P'jl<2)-U*pNl{l))/{ll-l.O) 

3 = ( c I ^> «» ( 0 + I ) ) • A ( N - 2 ) 

) G '< 1 ) = D n I 1 > * C 3 « { U 1 . 0 ' • B M 0 ( 3 ) 

;i ;? I ? ) = 0 A ( 4 ) • c 3 • P M 1 I 3 ) 

POO I 1 I = p:j0 I 2 ) 

P0in2)=P'j0‘3) 

P‘|4 i I ) sp-ji < 2 ) 
p *1 * 2 ) s p'l I < 3 ) 
in cnOTl'JUF. 

0G‘ 1 ) = jr, ( 1 » /t- 
.1-5 1 2 ) =0;i< 2 ) /r; 

[;GI3)=G,0 

- 0 
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Sl'MHO f ’Ifi ’ir r. ‘M X » OG I 'If ) 

1 I MEM i I N < » ijG ( i ) , Sr. ( 3 ) 

■) ! M5‘I3 : O'l 3 ( '’3 » ''<• I * SmL < ?3 I » GMl.( ?3 > 
COMMO'l/bGOtinn/AMM* 273 ) 273 ) 

COMMOU/fi AD/P. 1 * R 
COMMO!l/GRE/Rr 
S r I = X f iU / P, 
c F 1 aq i/f: 

!|Ml 

$')!,( 1 > B 0 , l» 

sml < 2 ) = x { 2 ) /n 1 

CML < I ) nl . 0 

c,u. ‘ 2 )*- 5 XU >/r. 1 . 

1)0 I Ib3,NMI 

SMI. « I ) = SMU < I - I ) *CML I 2 > *CML < I - 1 ).»SML < 2 ) 

1 CML < I > = -SML I t - I ) »SML ( 2 ) ♦OML* fl ) *CMU 2 > 
UM*tl,3 

i);i = » . 0 
R( i » I )= 1. 0 
P < 2 • L ) s s F ! 

P ( 2.7 2 ) = C F I 
tlM 2 s-lMl + l 
00 ?. lr: 2 »NM 2 

2 p ( I » T ) “*0 . 0 

)0 ■» I = . 3 »NM 2 

3 f> ( 2 ; I ) r 0 , 0 
1)0 4 MB.'., MM2 
i)n-M- 1 

0 0 M = 3 ,MMi 

.)= 1 ■?. . 0»UM“i . 0 ) 

IF.<M-M) 5 , 6 , 7 
5 P(MfM>“0.0 
Go To •'* 

o P ( M » M ) BU#Cf- ,t 41 p ( f,'-l » M - U 

i TtO_Tq._/, 


T }S f ^4 • t ' i ’.t ( I ; ! * r ' • . 

u CONTI lUE 
■ IUl')H = Rt;/n 
l)fl =0 . 0 

00 » '-l.i 

3 0G( I ) 

00 9 M=3,HMI 
UN- M 

no JLO !=l».s3 
10 s r. n ) s 0 . 0 
IJM*0 . 0 
1)0 1 I 11 = I » N 


•MArt=M f 

1 F t M 3 . L T . 1 ) . * 

Gf) TC 

30 12 JS=1*NS 
12 MAObmAC- ^JS *2 


13 


,11 


1 3 


4 


" V I'" U I ' 0 ' fl r- j tli C ! ’ • ) » / < U ” U • 
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C 1 = AMM.1 MAU ) ( M ) ♦PmMC nAI^ ) •SMj. ( M > 

um = m 

MS I U't- ; . ) *P ( tl , M ) 

S;'. ( 1 ) = sr» ( 1 ) *0 l«p { N , M> 

.SG‘ 2)=r,C'( <l >*Cl*<P<M,M*l)-f*<x‘3)/Rl ) ) 

Sr, t •' ) 3 i(, I • ; ; - A MM <IJ AB > • SML < M) *&MM < N Aft ) *CML ( M ) )*H 

02 = Kn0P, •••’■) 

;) 0 ‘ ). ) = D G I 1 > ♦ c 2 ♦ U tJ • S G ( 1 > 

■')/; «2)=DG(2)*c2*SG<2) 

00 ; 3 ) 5,00 ( 3 ) ♦gP^SG < 3 ) - 
OG I 1 ) = Dfi ( 1 ) /P. 

0 G < ?. 5 = 0 R ( ?. » / R 
0 0 ( A ) = 0 f, ( .3 ) / p 1 
H t. T U R .M 

K M I) 
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sii»f^nu7 inF. orr. ( X » Ofi ) 
x<3».rir, (3) 
Cl.MMnN/RAl»/R • , R 

cnM;\.v)/CR^./Rn 

‘*<nMM0N/r.Aa,'/ft;>?,h22, a 30, A<*0 
ci = Hf-/n 
n 2 » 0 I • C 1 
C 3 " 0 ,•» • c 1 / R 
c^. = r.3*n.'./R 
ui = )(t 3 )/r 
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1) 3 = X ' 1 • X < i > - y t ? >■ * X < ;> ) 

m? = hm 

R ix=m^f. 1 
tl<‘ = U3/|', ,\2 

= 2>/pi2 

jr. < l 0«X( 1 )*x( 2)«»P22/R2'»Cl*Ul 

• * ' 1 n , 0 i. !l 2 - , !J ) 4, A 7 r j. ( 1 7 5 , 0 * IJ 2 , >j 2 - 1 3 0 . 0 • u 2 •*■ 1 5 t 0 » • C2 • A <► U X S 

n r, t ?. ) c n «. « ( - 0 . 0 I ' 1 ^ 4 « J J 2 ^ 2- ^ g ♦ , , 5 , e., 2 2 J * J /5 « Q 1 

• •<‘^.0*u2-l,0)4.;;Ao.r;,5*c2*Ul*'<7tO«‘U'2*3,0)«AAO) 

Dr. o ) sc , 0 * / -r , ii * 11 ! “ A 2 2 + U4.62 2 ) 

K E T U K iJ 
ENP 
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List of Subprograms by Marne. 


Name 

Index 

Page * 

A8S£L 

fins 

?9 

ADEN 

DO^ 

50 

1 AGIGA 

BOP 

?5 

AGiOAO 

Br\f. . 

?5 

AMBAR 

■ D03- 

50 

CONOR 

Anr> 

T5 

CONST 

Ar)j 

TT 

1 datday' 

BOT 


DAYDAT 

sn? 

20 

DEQ2 

Enj 

«^T 

I DEQ?> 

EOT 

5T 

1 DEO A 

£ni 

5T 

j DENS 

DO’> 

45 

! SLABS 

Qnn 

P'7 


p o ^ C t & a/ T c- *: 7 " pw c e 
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4 ^- 


NAME 

INDEX 

PAGE 

gcltLn 

5 TO 

37 

GEOOaC 

B TO 

?3 

GRAY 

DO? 

SO 

HEIGHT 

Bin 

S3 

MMSSEC 

8 0^ 

?? 

'RO ■ ' 

DOT 

A? 

&EGHMd 

BOA 

?? 

SELENA 

COT 


STl' •• 

B0*^- 

?3 

’ STTIME 

Bnu 

?3 

SUN 

CO? 


TLOCAL 

DO^ 

SO 

TRDATO 

5 09 

30 

VHAiA 

i 

D 

A3 
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